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Abstract 



i-C ' Considering the possibility to build an e^e collider at the energies around Z-boson resonance 

o ■ with a luminosity so high as £ oc lO'^^cm^^s"^ (even higher) and the abilities of a modern 

^S^ , synthesis detector, we systematically calculate the exclusive two body processes: e^e annihilates 



into a heavy quarkonium and a photon (initiate state radiation i.e. ISR is involved), at the 
energies around the .Z^-boson resonance. Since the couplings of Z-boson to quarks contain axial 



> 

cn 

t^^ . vector as well as vector, so here the produced heavy quarkonium may stand for a charmonium 

^, 

j_^ , such as J/ip,'r]c,hc,Xcj ■ ■ ■ and a bottomonium such as T,r]f),hb,Xbj ■ ■ ■ respectively. If we call 

f— ^ ' such a collider with so high luminosity and running around the Z-boson resonance as a Z-factory, 

then our results obtained here indicate that experimental studies at a Z-factory about the various 



k> I heavy quarkona (their ground and excited states) via the two-body processes, especially, the 

production 
advantages. 



C^ ' production of the bound states with quantum number J^*^ = 1 via ISR, have outstanding 
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The production of a heavy quarkonium, such as a charinonium or a bottomonium, at high 
energies not only may be used to test the high-energy behavior of quantum chromodynamics 
(QCD) as well as the interplay of perturbative and non-perturbative phenomena in QCD, 
but also may be used to search for suitable sources to produce experimental example for 
studying the properties of the heavy quarkonia. Comparing with the hadronic colliders such 
as TEVATRON and LHC, an e^e~ collider has many advantages except the low rate of the 
production in studying the production and/or in using the produced heavy quarkonium as 
an experimental study sample. Recently the technique progress on e~^e~ colliders motivated 
by ILC has achieved, that the luminosity may be raised so high as £ oc 10^^cm~^s~^ (even 
higher). Moreover, if such an e~^e~ collider further runs at Z-boson peak energy, the res- 
onance effects at the peak may raise the production rate to several magnitude order high, 
thus the shortcoming (the low production rate) can be overcome much for the mentioned 
purposes. For convenience, we will call as a collider as a 'Z-factory' later on^. As estimated 
in Ref.[l|, in terms of GigaZ, one can perform experiments on the basis of an example more 
than 10^~^'' Z-events, that enlarges the number of Z-events by 3 ~ 4 orders of magnitude 
than that collected in LEP-I. Thus we think Z-factory will be able to open new opportuni- 
ties not only for high precision physics in the electro-weak sector, but also for the hadron 
physics including that of the heavy quarkonia. 

Considering the abilities of a modern synthesis detector for photons, and the fact that 
the energy-momentum is monochromatic in a two-body process, who's final state contains 
two body only and incoming energy is fixed, in the present paper we would like to focus the 
two body processes for the heavy quarkonium production as below: 

e^{pi) + e-ip2)^liP3) + HQQiP) (1) 

where the heavy quarkonium state Hqq stands for J/ip, rjc, he, Xcj {J = 0, 1, 2) for charmo- 
nium and stands for T, rjt,, hf,, Xbj {J = 0,1,2) for bottomonium respectively^, especially 
at the energy of Z-peak. Since the produced heavy quarkonium in the two-body processes 
is monochromatic in energy-momentum, so experimentally the heavy quarkonium may be 



^ In the literature (Ref.[l| and references therein) running the prospective high-energy e~^e colhder is called 

as GigaZ. 
^ In fact, here Hqq can stand for higher excited states too, but in the present paper we restrict ourselves 

to consider the low-lying excited states only. 
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FIG. 1: Typical Feynman diagrams for the process e^{pi) + e~{p2) — ?• 'y* /Z^ — )• 7(^3) + Hqq{P), 
where Hqq equals to J/ip, rjc, XcJ (J = 0,1,2) for charmonium and T, rjh, XbJ (J = 0,1,2) for 
bottonionium respectively. 

comparatively easy to be identified, as long as the investigated event-sample is big enough. 
Thus how great the cross sections are crucial for experimental studies of the heavy quarko- 
nia. Moreover to observe the production of the heavy quarkonia at Z-peak via the two body 
process, the resonance effects i.e. the contributions from 'Z-boson exchange' become greater 
than those from 'virtual-photon exchange' owing to the Z-boson propagator comes to its 
mass-shell, and 'Z-boson exchange' is of weak interaction (without parity and charge-parity 
conservation), hence various states of quarkonia with different parity and charge-parity can 
be produced quite equally. To the leading order calculation, there are four typical Feynman 
diagrams FIG. la, FIG. lb, FIG.lc and FIG. Id for the processes (since each typical Feynman 
diagram in FIG.l contains two diagrams: a virtual-photon exchange one and a virtual Z- 
boson exchange one, so the total number of the independent Feynman diagrams represented 
in FIG.l is eight in fact) if the the heavy quarkonia in the final state have quantum numbers 
as J^'" = 1 , whereas there are two FIG.lc and FIG. Id only (for the same reason as in 
the case for the heavy quarkonia J^'" = 1 , the number of independent Feynman diagrams 
are four) if the the heavy quarkonia in the final state have quantum numbers different from 
J^*^ = 1 . As matter of fact, FIG. la and FIG. lb are those of ISR (here one may see that 
ISR produce the states with quantum numbers J^"" = 1 only). Note here that due to the 
fact that a heavy-quark propagator via two vertices connects with the heavy quarkonium as 
shown in the diagrams FIG.lc and FIG. Id, but not as that shown in the diagrams FIG. la 
and FIG. lb, a photon directly connects to the heavy quarkonium via a single vertex, then 
the heavy quarkonium, involving orbital angle- momentum excitation (L > 0), can be also 
produced via this mechanism depicted by the diagrams FIG.lc and FIG. Id. 

The calculations on the processes, no matter for the ground state or for the excited states, 



need to establish the exact relations between each term of the amplitude (each Feynman 
diagram) and the nonperturbative matrix elements in NRQCD framework|2], or equivalently, 
the relevan Bethe-Salpeter (B.S.) wave function (derivative of the wave function) at origin 
for the 5'-(P-)wave production respectively in B.S framework 3|. 

Generally, according to Feynman rules, the term of the amplitude relating to the con- 
cerned Feynman diagram, can be written as 

M, = Cv,ip,)T,kUeiP2)t J -^,Tr[x{P,q)T2k] , (2) 

where x(-P) q) is the B.S wave function for the concerned quarkonium, C stands for an overall 
parameter, which is the same for all the Feynman diagrams. Fi^ stands for the structure of 
the fc-th Feynman diagram, and is for the electric line, which includes the relative string of 
Dirac 7 matrices and the corresponding scalar part of the propagators, r2k stands for the 
similar structure of the heavy quark line. For convenience, the intermediate photon or Z 
propagator is put into Fi^. 

As for a quarkonium in S'-wave state, in the non-relativistic approximation, x{Py 1) can 
be written as 

X{P. q) = 0(g)^(a75 + /3e(^.))(P + M) , (3) 



where P, M c:^ niq + niQ are the momentum, the mass of the heavy quarkonium respectively 
and q is the relative momentum between the two constituent quarks. In this paper, we use 
a to denote the contraction between the Dirac 7 matrix and a momentum or polarization 
vector a, i.e. we use a instead of (i. a = 1, /3 = for the pseudoscalar ([^^o]) quarkonium 
and a = 0, /3 = 1 for the vector (['^5'i]) one. The wave function 0(q'), the radial part 
of the momentum space, can be related to the space-time wave function at origin by the 
integration, i j 72^0(?) = V'(O), where |V'(0)P = |-R(0)p/(47r). Then the 5'-wave amplitude 
can be simplified as 

1 



Mk = Cve{pi)TikUe{p2)^{Q)Tr 







-^(«75 + /3e(^.))(P + M)F2, 



(4) 



where F^^ is defined by the following Eq.Q. 

As for a quarkonium in P-wave state, in the non-relativistic approximation, the B.S. wave 
function x(-P) q) = x' P.jJz](^p^ g^ ^g^^ l^g written as 



^[<--'P...l(P,,)^ E p^ . v,(,) . (e^ . ,) . 
5., A, A' 4mQmQ 

- "^q) ■ ('55,o<5s.,o75 + 5s,i5s,,\>e^') ■ (qg + mq) ■ (lA; SS^\JJ^), (5) 
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where qq and qq are the momenta of Q and Q quarks inside the quarkonium: 

iriQ ^ _ mq 



qq = aiP -q, qQ = a2P + q, ai 
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(6) 



mq + mg mq + mq 

For convenience, we introduce g^ = q^ — qi!', qit = ^^j-^P^ and q\\ = \qit\, hence accordingly we 
have d'^q = dq^^d^q±. \I/(g) stands for the 'P-wave scalar wave function', and (lA; SS^lJJ^) 
is the Clebsch-Gordon coefficient for L-S coupling. The spin structure of the B.S. wave 
function x^ ^'^'"'^'(g) defined in Eq.(l5]) is of the lowest order in q) for the P-wave state 
production. 

The derivative of the wave function at origin V'qIO) i^ coordinate representation relates to 
the P-wave scalar function "^{q) in Eq.Q under the so-called instantaneous approximation 
by the integration: 



dq\\d^qj 



■'°*(')(''-')='/Sf'°H"sj^ 



(e^.g^) = ie-^"^'(0) 



(7) 



g± 



„2 ' 



and 
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■jj^) stands for 



where |'?/''(0)p = |P'(0)p/(47r). Here qj_ is the unit vector 

the P-wave 'scalar wave function' in the sense of potential model. Substituting Eq.([6]) into 

Eq.dS]), we obtain 



'{v7JI^~^ + ^^ ■ ('^^'O'^^-o^s + Ss,iSs^,x'e^') 



-( 
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+ai{-P + M) • {5sfl5s-^fll5 + 5s,i5s^,x>e^')qi\ + C(gi)} 
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where 
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P + M)- {6s,oSx',ol5 + Ss,iSs^,x'e^') , 
a2Y{-P + M) ■ {6s,oSx',ol5 + Ssa^s^^x'^^'] 



\4mqmQ^ 

+ai{-P + M) ■ i6s,o6x',ol5 + 5s,i6s^,x'e^'h'' 



TABLE I: Functions C -^ , EJi '^ ^ and F^ Pj\J^ for a particular P-wave heavy quarko- 

nium states, where e^ and ej^l, are polarization vector and tensor respectively. 
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To leading order of the relativistic approximation and according to Eq.([2]), the next step is 
to do the expansion of the r2k about q^ up to 0{q^) for the P-wave production, i.e. 



^2k = ^k + r^, • q, + 0{q') 



(9) 



where T^^ and V2k are independent of q^. Substituting the above equations into Eq.([2]) and 
carrying out the integration over d!^q = dq\\(Pq_[_ and with the help of Eq.([7]), we obtain 



Mu = Cve{p,)r,kUe{p2) E5.,A,A'(1A; S S^\J J^)^' {0)e'^ ■ Tr [A'ss. " T^fc + ^5^1 " ^2k\ ■ (10) 

For a specific P-wave state, summing over the explicit A and A' for the Clebsch-Gordon 
coefficients, Eq.lTTOl) can be further simplified as 



M, 



SJJ:, 



CvM)ri,u^iP2)c''''''^i^'mr fi^r^^^^i^^ ■ r^, + pP-^^^]^^ . Tl 



:ii) 



where the overall factor C ^-^ and the functions E[ ^-^'^^ and F^ ^-^^^^ are shown in 
TAB HI 

To do the computation we adopt the 'linear polarization' for the heavy quarkonium vector 
and tensor states accordingly, and their explicit formulation can be found in Ref. [4|. The 
amplitude of all mentioned processes can be easily generated by using FDC^. To simplify 
the amplitude, firstly we establish a complete set of ' basic spinor lines', such as that they 
are constructed by the multiplication of Dirac 7 matrixes in a certain way as P1P2P3 ■ ■ ■, 
then we consider them as bases to expand every terms of the amplitude and sum up all the 



The relevant Fortran program for the mentioned processes is available upon request. 
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FIG. 2: (color online) Total cross sections for the processes e~ +6"*" — )• j + Hqq versus the collision 
energy. The red solid, the black dotted, the blue up-solid-triangle, the green dash-dotted, the red 
dashed and the down-hollow-triangle lines stand for QQ in ^Si, ^Sq, ^Pq, ^Pi, ^P2, ^Pi respectively. 
The left figure is for charmonium and the right one is for bottomonium. 

terms according to the expansion (the coefficients of the 'bases' are summed respectively). 
By this way, the amplitude and its square become quite condense, thus the efficiency for 
computing the amplitude squared is raised greatly. This kind of treatments can be carry 
out quite automatically in FDC package nowjsl]. 
When doing numerical calculations, we adopt 



mz = 91.187Ge1/ , rrib = A.QGeV , rUc = 1.5GeV , 



(12) 



and 



?(cc) I 

hs 



\R\''.^KO)\' 



O.SWGeV^ , |i?fsHO)r = QA47GeV^ , 
|i?g)(0)|2 = 0.075GeV^ , |4p (0)1^ = lAUGeV' (13) 

for the parameters of heavy quarkonia [6]. In order to have some idea about more samples 
to produce a 'particle' with such quantum numbers as J^^ = 1 via the ISR in addition 
to the ground states J/ip and T, 

To present the results, the dependence of the total cross sections for the processes e"*" -|- 
e~ — )■ 7 -|- Hqq with Hqq for charmonium and bottomonium on its collision energy is 
drawn in FIGEl cind the dependence of the differential cross sections for the processes 
e~ -|- e"*" — )■ 7 -|- Hqq on cos a, where a stands for the angle between the quarkonium and 
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FIG. 3: (color online) Differential cross sections for the processes e~ +6"*" — )• 'j + Hqq versus cos a 
at a C.M.S. energy as Z-mass. The red solid, the black dotted, the blue up-solid-triangle, the 
green dash-dotted, the red dashed and the blue down-hollow-triangle lines stand for QQ in ^Si, 
^5*0, ^Pq; ^-Pi) ^^2) ^Pi respectively. The left figure is for charmonium (the dotted line and the 
blue down-hollow-triangle almost emerge together almost) and the right one is for bottomonium 
(the red dashed line, the green dash-dotted line and the blue down-hollow-triangle emerge together 
almost). 

TABLE II: Total cross sections for the different quarkonium states at the collision energy y^ = mz- 
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the beam, is drawn in FIGjSl From FIGj2l one can see clearly that as expected there is 
a peak at ^/s = mz for the quarkonium states of ^^o, 'Pi and ^Pj etc. It is because 
that in these processes the s-channel diagrams of Z-boson exchange FIG.lc and FIG.ld 
become dominant, when the Z-propagator in the diagrams approaches to mass-shell. We 
should emphasize that the Z-boson coupling is not only of vector, but also of axial vector, 
so the states with various quantum numbers, as long as they are lighter than mz, may be 
produced in terms of the s-channel diagrams of Z-boson exchange, and there must be a peak 
at y/s = mz- In contrary, if the production is not at Z-mass energy but at an energy much 



low, then of the s-channel diagrams FIG.lc and FIG. Id, Z-boson exchange ones become 
ignorable and only the ones of 7 exchange play a role, so that there is a vector coupling 
only, and quite a lot heavy quarkonium states without suitable quantum numbers cannot 
be produced in terms of the 's-channel mechanism'. 

However, for a hidden flavored 'comparatively light' quarkonium {TriH ^ ^^z), when 
the quarkonium is a ^Si or '^Di state i.e. that with quantum numbers as J^*" = 1 , then 
the t-channel photon-exchange diagrams, as shown in FIG. la and FIG. lb, will be dominant 
over the other diagrams (including FIG.lc and FIG. Id even the Z-propagator being on-shell 
and the t-channel Z-boson-exchange ones in FIG. la and FIG. lb). In fact, the 'dominance' is 
also the nature of ISR essentially. Thus the production of a J^*" = 1 heavy quarkonium 
is very large and without such a peak behavior around the energy of Z-boson as emerged in 
the production of some other states. To show the fact, the dependence of the cross-sections 
for producing this kind of quarkonia is also drawn in FIGJ21 

Moreover, to make comparison, we put the total cross-sections for the production of the 
all low-lying states at Z-factory in TAB|TT1 We would like to note here that if varying 
a/s = m^ ±0.5 GeV slightly, the results will bring ±1% difference in the total cross sections. 

From FIG 131 one may see that the distributions for the production of the heavy quarkonia 
are quite fiat in directions, although in the forward and backward directions they arise in 
certain aspect, thus there is not essential loss on observing the heavy quarkonia, if putting 
a necessary cut around the beam directions. 

To see the role quantitatively that ISR may play in the two-body final state processes 
at Z-factory, we assume a 'X(4260)-like' particle with quantum numbers J^*" = 1 and 
calculate its production precisely. For convenience, we denote the 'X(4260)-like' particle as 
X^5 below, and for definiteness we further assume its mass is ^tix^^ — 15 GeV (^^x^^ > '^t), 
its decay constant is one third of that of T: /x^- = |/t and its main decay mode is similar to 
that of the particle X(4260) [TJ, i.e., X^g — )■ Tvrvr, but not X^,i — )■ J jipixi^ so same as X(4260). 
Finally we obtain the total cross-section of producing X^i at Zq peak is 1.898 ■ 10"^p6; while 
as a comparison, the total cross-cross section of producing T is 5.650 ■ 10~^p6. 

Summary: From the results obtained here, one may see that the production of a particle 
with Quantum numbers J^'^ = 1 via the two-body final state processes e~^+e~ — )■ '-f+NgQ 
concerned here is of ISR essentially and the total cross sections for the production of the 
heavy quarkonium states at Z-boson peak via the two-body final state processes can be so 



great cr ^ 1 ■ 10~^pb in magnitude, when Hqq is a heavy quarkonium with quantum numbers 
J^*^ = 1 ; so great o" ^ 1 ■ 10^'^pb, when Hqq is a 'X(4260)-hke' particle with quantum 
numbers J^*" = 1 as assumed in the above; whereas the total cross sections at Z-boson 
peak for producing the quarkonium states Hqq with the other quantum numbers (not J^'-' = 
1 ) via the two-body final state processes (without ISR) can only be a ~ 10^'^ ~ 10~'^pb in 
magnitude. Therefore, one can conclude that numerous events of Hqq with J^'-' = 1 (the 
heavy quarkonia and even a 'X(4260)-like particle such as X^g as well) may be produced at a 
Z-factory. Also at such a Z-factory, hundreds, even thousands events of the heavy quarkonia 
Hqq with different quantum numbers from J^'^ = 1 may be produced. Considering the 
advantage of the mono-energy characteristic for the produced photon and heavy quarkonium 
both in the two-body final state processes, as well as the clean environment in e~^e~ collision, 
it seems that it is worth more efforts quantitatively to estimate the precise advantages in 
studying the production and the spectrum of the heavy quarkonium systems additionally 
with specific condition of a synthesis detectors at Z-factory. One can also conclude that it 
will be better if the luminosity of the Z-factory may reach to so high C > 10^^~^^cm~^s~^ 



.0 



for the heavy quarkonium studies. Moreover in Ref.J8|, we report additional advantages in 
studying the heavy quarkonium properties at a Z-factory, but via some processes else. 
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